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PREFACE 


The large variety of ablation markings obstirved on recovered 
tektites has lead to the swarm wake model first proposed in the 
preceding investigation. Herein, further considerations are 
presented in support of this model. Quantitative assessments 
indicate that wake shielding might indeed have provided for sub- 
stantially less heating than would have been experienced by a tektite 
entering an undisturbed atmosphere along a similar trajectory. 

For the case of strong wake shielding it is even possible that 
the surface temperature of a falling tektite had bcrely reached 
its melting point. It is argued that in the distribution of 

tektites there is a size band (near R = 0.5 cm) which is least 

* 

susceptible to melting. 
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ATMOSPHERE ENTRY OF A SWARM OF TEKTITES 


1.0 INTRODUCTION 

In a preceding study^^^ it was stated that several re- 
covered tektites exhibit surface features which indicate that 
these tektites had sustained very little ablation during their 
passage through the earth's atmosphere. This observation is 
rather puzzling since various calculations concerning the moat 
likely tektite trajectories imply heating rates compatible with 
substantial ablation. A reconciliation was proposed by means of 
the swarm hypothesis# which states that the lead peripheral 
tektites bore the brunt of aerodynamic heating upon entry# and 
that the bulk of tektites in the wake enjoyed partial shielding 
at the expense of the leaders. This shielding was embodied in 
a shock envelope model which implied mitigated entry conditions for 
the trailing tektites through a decrease of effective entry 
velocity and appropriate f;' \.?iges of the ambient density. 

The objective of the present investigation is to analyze the 
swarm hypothesis in greater depth in order to assess its viability 
as a solution. It is concluded that such a circumstance could 
indeed have reduced tektite ablation to the degree actually 
observed . 
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2.0 THE MOTION OF A TEKTITE IN THE SWARM WAKE 


The ambient conditions of a tektite entering in the wake of the 
hypothesized shock envelope would have been composed of the effects 
of a sequence of preceding tektite wakes (randomly spaced) such 
that each successive wake would have amplified the pertinent 
effects. In this light it is difficult to assess quantitatively 
the decrease of tektite ablation. However, some appreciation of 
the phenomenon should be attained through consideration of the 
simpler case of single w^e shielding. Therefore, we initially 
restrict our attention to the descent of two tektites, one of 
which follows the other somewhere in its wake. 

Experimental evidence of hypersonic wake characteristics 
has appeared mainly in the form of Schlieren and holographic 
photography and recently in the measurement of wake velocity 
profiles by means of a sequential spark technique. Beyond 
the near wake region the velocity defect distribution attains 
a similarity form that may be approximated by a Gaussian representa- 
tion such that the centerline decay approaches a power law of 

(2)-(5) 

exponent -1 in approximation . Evidently,' there rs no 

direct data available concerning the density structure of a 
hypersonic wake; consequently, we resort to the following 
theoretical argiament to obtain the model required for subsequent 
analysis. 
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2.1 Average Momentum Deficit in a Hypersonic Wake 


Consider a tektite descending vertically in the earth's 
atmosphere. During the period in which the gravitational force 
is negligible, conservation of momentum may be expressed as: 

mvlt) = mv:, - U) 

vet) 

where the integral term represents the momentum gained by the 
air, and the voluitie of integration includes only the air dis- 
turbed by the interaction. The above may be re-written in the 
simple form: 

'Vl*} _ I _ Vt-t) 

mvio ( 2 ) 

where 

V 

This integral represents the average momentum density gained by 
the disturbed atmosphere. In a hypersonic entry the disturbed 
atmosphere is contained in a vol\ame bounded by the conical shock, 
so that; 

Vtt) = ^irta^CB) lUi) 

3 (3) 

where represents the total distance penetrated by the 

o 

tektite, and ©• is the cone half -angle, which is typically 12 . 


Combining (2) and (3) « one obtains: 






Xio0d(&)r 

•9 


[ 



( 4 ) 


Equation (4) provides an expression for the average momentum 

deficit in the wake as a function of the distance downstream of 

the tektite. The velocity of the tektite as a function of 

(6) 

altitude is given approximately by: 



<5) 


where the atmospheric density is assvimed to decrease exponentially 


with altitude: 






an^ 

<»<- = f Co ^ 2.09 X 10”'^ 

8 pi 

~ = 150 for R = 1 cm 

pR. 

The velocity profile for such a typical entry is shown in 
Figure 1. Combining (4) and (5) we obtain an expression for the 
average momentum deficit of the total wake as a function of 
tektite penetration into the atmosphere: 


i*. 


ir 




TT 


tutLe) ( 




( 6 ) 
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Here# y = 0 denotes that earth's surface and (4 0) denotes 
the tektite entry altitude at velocity 'V'^ • 

We now use Equation (6) in approximation in order to obtain 
an expression for the average momentum decay downstream of the 
tektite at any point along its trajectory. Let s " y)<0 

such that b/yqq « 1. Then expand (6) to obtain: 



m. t Co 


This expression indicates that the average momentum deficit 

decays inver?: with the square of the distance downstream of 

the tektite. In combination with the previously mentioned experi- 
0 

mental observation that the velocity deficit decays inversely 
with distance# the following model is proposed to describe the 
average density and velocity decays in the wake of a tektite: 


= f’-[' ~ S+S.3 

sfe] 


S < Cf 

(7) 


Such that the momemtum deficit is in accord with thri preceding 


argument: 
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This model shall be used to describe the ambient conditions 


of a tektite following at distance s in the wake of a preceding 
tektite. The model is not accurate in the near wake (s — *>o) 
where the wake decay alters its character, but it appears to be 
adequate for present purposes. 


2.2 The Motion of a Shielded Tektite 

Consider two identical tektites following vertically to- 
wards the earth's surface, one following in the wake of the other. 
The objective of this section is to describe the relative motion 
of these two tektites during their descent. 


Even in the absence of a wake effect the relative spacing 
of the 44' j tektites will decrease due to the lagged decelera- 
tion of the second one. A relation for this separation is obtained 
by noting that since both tektites may travel identical trajec- 
tories (albeit time lagged) , the time lag of their passage at 
any point in the trajectory is independent of position along the 
trajectory. Prom this fact one obtains: 


5 ^ 


( 8 ) 


where tT is the value of tektite velocity at some position of 


the trajectory between the two tektite locations at any given 
altitude. As the initial tektite velocity of 11 km/sec decreases. 


the separation decreases commensurately . Therefore, the initial 
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tektite separation may be larger than the extent of effective 
wake influence, but the second tektite may close into this domain 
at some later time during the trajectory. 


When the motion of the second tektite is influenced by the 
wake of the first, the following equations may be used to 
approximate their respective motions: 


It 


2t 


( 9 ) 





( 10 ) 


where 


« I t 


C a' 

0 


For the second tektite the 


8 7, « 

ambient density and velocity are represented as being perturbed 


from the undisturbed atmosphere by the amounts & and S 
respectively. We define the relative separation and velocity 
between the two tektites as follows: 

s ® yj - <0. 


w s ir . TT 
2 1 


Subtraction of Equations (9) and (10) leads to the equation of 
relative motion between the tektites: 


iJt 

dLt 


-- cse ( ^ -#-£)( u> “ •' — -f (u) 


( 11 ) 


where 



-0-rj ± 


t 


f - E — 



The previously mentioned wake model suggests that f-f; = X 
Since and UJ 2 functions of s (and therefore of t) , 

the solution to (11) is complicated. However, a qualitative 
understanding of the solution is readily apparent. For an initial 
relative velocity lying between the values and UJ ^ the 
forcing function f(u>) is positive. Therefore, the relative 
velocity will increase until it reaches the value ^ 2 ' whereupon 
further increase of u) is not possible. The second tektite 
will thus experience relative acceleration towards the lead 
tektite until it reaches the near wake with a positive relative 
velocity. Thereafter its momentum will likely carry it past the 
first which will then ta^'e the role of the shielded tektite. 

The analysis also demonstrates that even if the relative 
velocity were initially negative, ^e second tektite could 
catch up to the first. Therefore, net only is there a tendency 
for trailing tektites to catch the leaders, but there is also 
the tendency for all lead tektites to remain clustered near the 
front. In this manner the peripheral tektites form a more 
densely populated shield which sequentially creates a larger 
and more effective wake for stibsequent entries . 
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It is important to provide an estimate of the time required 
for a trailing tektite to catch the leader. Since the relative 
velocity between the two tektites is positive even without the 
wake effect# let us assume it takes the value 2 initially. 
From equation (11) the solution is simply: 

a; S 

Then 

S = u)gt + Si 

t - 


The value here depends on the strength chosen for the wake. 
Approximations are now based on the wake model pre <=:nted earlier. 
According to Equation (7) let: 

£ sr r = - C 

( 2 ) 

Experimental data for the velocity decay in a far wake in- 
dicate the values: 


c 

ZR 


20 


2R 


= 200 


Therefore, for a tektite traveling in range 200 4 S/2R < 1000, 
one computes the average value: 

£ a 20 - _ __ _ 1 

^ “ 200 +600 “ 40 
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Conclude: 




Therefore, an initial tektite separation of 10 m. will have 
closed in a descent of approximately 1 km. for the pair. 

2 . 3 The Effect of owarm Entry into the Atmosphere 

In a swarm of tektites entering the atmosphere the walte 
effects are further complicated by random distributions of tektite 
size, location and relative velocity. Relative decelerations 
are strongly dependent on tektite sizes and geometries. However, 
a tendency has been demonstrated here that lead tektites will 
bunch near the swarm perimeter. Tektites entering in the ensuing 
wake will experience heating less than that of a normal entry 
into the unperturbed atmosphere. 

Tektite showers were apparently events of considerable 
magnitude since some strewn fields are a few thousand miles long 
and contain a great deal of mass. The initial phases of such 
an event probably would have put the atmosphere into violent 
motion such that the later phases may have occurred with ambient 
atmospheric conditions greatly perturbed from normal. The wake 
flow accompanying the descent of the lead tektites would have 
compressed the atmosphere immediately behind the shock envelope 
and left the higher atmosphere relatively depleted of air. 
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The ensuing high and low pressure areas might have caused the 
generation of a toroidal vortex in the atn^sphere similar to the 
kind observed in a liquid subsequent to the impact of a droplet. 
Although the velocity perturbation of such a vortex would 
likely have been insignificant to subsequent tektite entries, 
the change in the atmospheric density profiles could have been 
considerable. The net effect of the motion would have been 
to reduce the height of the effective atmosphere while causing 
the lower altitudes to become denser. Such an alteration would 
have delayed the moment of peak heating. 


- 12 - 


3.0 THE HE^TIMG OF A TEKTITE IN THE SWARM WAKE 


(Aether a tektite enters in a wake or into an undisturbed 
atmosphere, it still must lose the same amount of kinetic energy. 
This energy is partitioned among several modes which are expressed 
by the energy balance: 

i W. ( •v’' ) sr ^ Cj CaT)^ +. ntC, (aT)^ ^ 

Several of these terms are negligible in the overall balance. 

(i) The final tektite kinetic energy is small compared 
with its initial value. 

(ii) Momentum considerations imply that the gain of 
kinetic energy by the atmosphere is small compared 
with the initial tektite kinetic energy. 

(iii) If a typical tektite were to absorb all of its 

kinetic energy internally, then its temperature would 
o 

rise by 50,000 K. Since some of the tektxtes do 

o 

not even reach the melting point (approx. 2,000 K) , 
tektite heating does not account for more than a 
couple percent of the original entirgy. 

(iv) After the completion of recombination of dissociated 
gas in the atmosphere, chemical energy levels return 

to their initial values, so there is little net change 


of this term. 



Therefore* the bulk of tektite kinetic energy goes into 
heating of the atmosphere* radiation loss* and to a lesser 
extent into the latent heats of melting and vaporization of 
whatever tektite material is ablated. Since the total heat 
absorbed by a barely ablated tektite is small in comparison to 
the other factors* it is possible to r^^duce this amount signifi- 
cantly while the major balance of energy is only slightly altered 
in proportion. Here we claim that entry in a wake may cause a 
significant decrease of tektite heating over that of an isolated 
entry* and that such a decrease would be absorbed by an imper- 
ceptible increase in atmospheric heating. 

3.1 Gas Cap Shielding 

One aspect of decreased tektite heating in a wake may be 
explained qualitatively as follows. Surface heating depends on the 
velocities of molecular impacts and on the gas density there 
since these quantities are uniquely related to the kinetic energies 
of impacts and their frequencies. In the rarified portion of 
entry the air molecules strike the tektite surface directly 
since the mean free path is large. As the density increases, 
the frequency of impacts increases, but eventually a gas cap 
forms in front of the tektite which absorbs much of the energy 
of further incoming molecules. The ensuing surface heat trans- 
fer is due to conduction in the gas for which the surface 
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collisions are far less energetic than the direct ones. Much 
of the energy absorbed by the gas cap is swept downstream with- 
out heating the tektite. Figure 2 taken from Reference 7 
shows heat transfer as a function of Reynolds nunber for a 
rarified gas flow. As applied to a tektite trajectory (at 
constant velocity say) the heating is seen to increase with in- 
creasing density until a peak is reached whereafter the heating 
decreases. The decrease is taken to be connected with gas cap 
shielding. The heating in a wake entry during the rarified 
regime is diminished over that of an isolated trajectory for 
two reasons. First, the relative velocity with respect to the 
gas is lower. Second, the duration of peak heating is shorter. 
In the early part of the trajectory the ambient density is lower 
than normal because the lead tektites have displaced the ambient 
air molecules. In the latter part of a shielded trajectory 
the ambient density is higher than normal because air molecules 
have accumulated in the preceding wakes. Over these con- 
ditions a shielded tektite is subjected to peak heating for a 
shorter time than otherwise. 

This mechanism also provides a means for the atmosphere to 
absorb a greater portion of energy than for the case of an 
isolated trajectory. Entry into a delayed but steeper density 
gradient results in a thicker gas cap. Although the volume 
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of the gas cap is larger, its temperattare need not be appreciably 
increased because dissociation and ionization processes limit 
the local temperature. A larger volume of hot gas would also 
result in more energy loss by radiation. 


3.2 Stagnation Point Heating 

Although the heating of a tektite during its descent through 

the atmosphere is a complicated phenomenon, a relatively simple 

expression appears to provide a good approximation for its 

stagnation point heating rate. The expression has been dis- 

( 6 ) 

cussed by Allan and Eggars , has been supported by laboratory 

( 8 ) 

data , and recently has compared favorably with data taken 

/g\ 

during an actual entry test program' ’ . In metric units this 
expression for laminar heating is: 


% 


1 / 

i^y^l ^ - \ Kui 

Ry** \foj ^ M fcm./;ec j CIR** 


( 12 ) 


where the tektite radius, R, is given in cm. 

Figure 3 shows the heating rate as a function of altitude 

for a typical tektite (R = 1 cm) entry into an undisturbed 

2 

exponential atmosphere. The maximvun value of 6.4 kw/cm 
occurs at an altitude of 50 km. It is easy to show that for such 
trajectories the value of maximum heating rate remains constant 
with a variation of tektite radius. However, the altitude 
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at which maximum heating occurs increases for smaller tektites 
according to the relation: 



From Figure 3 it is seen that the tektite experiences stagnation 

2 

point heating of more than 2 kw/cm for approximately 3.5 

seconds of its descent. The analysis of the preceding investiga- 

( 1 ) 

tion indicates that such a heating rate would raise the sur- 
face temperature past the melting rioint in less than 0.3 
seconds. Hence, a tektite entering the undisturbed atmosphere 
under these conditions would surely have melted superficially. 

In order to use the same expression in describing heating 
in the swarm wake, one needs to know how the density and relative 
velocity profiles as a fvmction of altitude had been altered 
from the undisturbed case. The ambience in the swam wake de- 
pends on three effects: (a) The number of preceding tektites 

to enter the atmosphere, (b) the separation downstream of the 
immediately preceding tektite and (c) the degree of motion of 
the whole atmosphere induced by the preceding tektites. Since 
the length of a typical strewn field j.s approximately 100 times 
the effective height of the atmosphere, one would expect that 
during the later phases of a swam e'W'ent the local atmosphere 
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would have been in violent motion. Of course, for shallow 

entry angles the wake effect would be complicated by the swarm 

geometry with respect to the atmosphere and the density gradient 

with altitude. However, the general wake effect during the 

moments of maKimum heating would have been as follows. The 

ambient density would have been higher and the tektite velocity 

relative to the air would have been lower than those during an 

undisturbed entry. In continuum portion of entry the heating 

would have been significantly reduced as seen from Equation (12) . 

For the sake of comparison with the undisturbed entry, suppose 

the local density had been increased by a factor of four while 

the local relative velocity had been decreased by a factor of two. 

Then from Equation (9) it is seen that the trajectory for this 

wake case would have been identical to the case of the undisturbed 

entry. However, the heating rate from Equation (12) would have 

been reduced by a factor of approximately 4.5 throughout the 

trajectory. Thus, in Figure 3 tlie stagnation point heating 

2 

rate would have exceeded 1 kw/cm for only approximately 1.75 

seconds of the trajectory, while reaching a peak of only 
2 

1.42 kw/cm . If we assume that the tektite (R = 1 cm) surface 

2 

experienced an average heat flux of 0.5 kw/cm and that the 
ten®)erature external to the tektite in the reduced Mach 
ntimber flow had been decreased to 4,000° K, then the previous 
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analysis (Figure 4) shows that the tehtlte surface would 
have reached 2«000^ K In approximately 1.5 seconds. We 
conclude that under circumstances similar to these a 
tehtlte may have entered the atmosphere at 11 km/sec and 
yet have experienced heating only sufficient to have brought 
the surface barely to the melting point during Its descent. 
Furthermore# since the wake has such a large Influence on the 
heating rate# It Is clear that a large variety In the amounts 
of ablation would have ensued among the members entering In a 
swarm. Indeed# this conclusion conforms with the observations 
made Initially In connection with the varle^ty of surface 
markings on recovered tektltes. 

3.3 Turbulent Heating 

The possible occurrence of tvirbulent flow may have 

further complicated the question of tektlte heating. However# 

( 10 ) . . 

Reynolds nuiriber considerations xndicate that the transx- 
tlon to turbulence would have been a marginal event. The Rey- 
nolds number peak In the trajectory occurs subsequent to the 
deceleration peak which occurs stabsequent to the laminar heating 
peak. Hence, if transition did manage to occur during the tra- 
jectory# it would have been at relatively low velocity. 

Thus, the decrease in heating rate after the laminar peak 
would have been more gradual if turbulence had been generated, 
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bat the peak heating should still have occurred during the laminar 
phase. However « the total amount of ablation could have been 
significantly greater for the trajectory involving turbulent 
heating because the tektite surface would already have beea 
molten as the turbulence commenced. The increased heating due 
to turbulence would have lengthened the segment of the tra- 
jectory for v;hich the surface experienced melt flow and 
vaporization. However, the swarm wake model complicates the 
question of turbulence inception because the Reynolds number 
depends on the product of local density and velocity. De- 
pending on the actual relative changes of these quantities, 
turbulence may have been either inhibited or enhanced. 


3 .4 The Dependence of Ablation on Tektite Size 

The results for entry into an undisturbed atmosphere 
indicate thst velocity, heating, and Reynolds number vary with 
altitude and tektite radius according to the following functional 
forms: 




= R 








( 14 ) 
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Thus, in this model, every tektite would pass through the same 
trajectory and heating profiles except that the profiles are 
shifted in altitude with the location of the peaks given by: 




It should be noted that tlie magnitude of the Reynolds number 
profile is amplified by the factor R , whereas the amplitudes 
in the velocity and heating profiles are independent of 
tektite radiu~ . 


In the light of these profiles we consider the rise in 
surface temperature for various tektite sizes. The short time 
temperature solution presented in Idie previous investigation 
is valid until the time at which the tektite center begins to 
rise in temperature. This characteristic time is naturally 


shorter for smaller tektites, and the condition is expressed by 






Assuming the center begins to be thermally involved when the 
parameter reaches the value 0.1, we obtain the relation: 




( 16 ) 
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Therefore# for tektites smaller than 0.5 on in radius, the 
center begins to rise in temperature less than 1 second after 
the heating is applied to the surface. A>fter this characteris- 
tic time is surpassed, the tektite surface temperature rises 
more quickly than the short time solution indicates because 
the center becomes a less effective heat sink in comparison 
to that of a larger tektite. The short time solution demon- 
strates that initially all tektite surface temperatures rise at 
the same rate (along each trajectory) provided the surface heat- 
ing rate is the same, which is the case for entry into an un- 
disturbed atmosphere. The surface temperatures of the smaller 
tektites rise faster than those of the larger ones only after 
their centers have become involved thermally. 


Since it takes approximately 5 seconds along each tra- 
jectory for the tektite to reach its point of maxim\am heating. 
Equation (16) shows that all tektites of radius greater than 1 cm 
should begin to melt at the same rislative point along td:ieir 
trajectories. Tektites smaller than 1 cm should begin melting 
progressively earlier. Since the remainder of the laminar 
heating profile is the same for all tektites, one wovald not ex- 
pect much difference in the ablation of the various larger 
tektites. On the other hand, the Reynolds number profile 
increases in magnitude as R . Thus, for increasing diameters. 
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transition to turbulence occurs earlier along the relative tra- 
jectories, and progressively more of the heating profile after 
the laminar heating peah is augmented by the turbulent rate. 

Under these circumstances one would expect the amount of 
ablation observed on various tehtites to increase with increasing 
radius. Indeed, Figure 6 of Reference 1 shows such a trend 
markedly . 

Furthermore, this graph tends to support the notion that 
most tektites have experienced wake shielding to some degree. 

An extrapolation of the data points suggests that many tektites 
of roughly 0.4 cm radius had not ablated at all (or had 
ablated very little) . If such tektites had entered the undis- 
tvirbed atmosphere at 11 km/see and h"d experienced only laminar 
heating, without the benefit of wake shielding, then the analysis 
had demonstrated that they should have ablated significantly. 
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4.0 CONCLUSIONS 




/ 


The analysis of tektite heating and motion has lead to the 
following conclusions: 

(1) No tektite should have remained superficially un- 
molten if it had entered the undisturbed atmosphere 
at 11 km/sec. 

(2) If a tektite had entered in the strong wake associated 
with the later phases of a swarm event, its heating 
could have been reduced to the point of avoiding sur- 
face melting. 

(3) The surface (and internal) temperatures of microtektites 
should have risen to the melting point more rapidly 
than those of the larger tektites. 

(4) The heating rates of the larger tektites should have 
been increased via turbulence, thereby making them 
more susceptible to ablation than tektites of 
intermediate sizes. 

( 5 ) Tektites of approximately 1/2 cm radius were least 
susceptible to melting because these optimized the 
trade-off between turbulent heating and the low thermal 
inertia of the smaller tektites. 


/ 

(6) Much of the previously puzzling scatter in ablation 
data may be attributed to the variability of heating 
in a swarm wake. 

In connection with the present conclusions it is interesting 
to mention a general trend observed in some strewn fields in 
Indo-China. Xn these distributions the tektites which exhibit 

( 11 ) 

most ablation are foimd near the peripheries of the strewn fields 
This observation is in agreement with the notion of wake shielding 
because tektites in the periphery of a swarm should have enjoyed 
less of the wake effect than those towards the center. 

Xn addition, the conclusions of this report raise implications 
connected with computed re-enactments of tektite trajectories 
in which ablation data is used to determine the trajectory. 

Since these computations assume entry into an undisturbed 
atmosphere, it would be most accurate to consider those speci- 
men which exhibit the most ablation. The fianged australites 

( 12 ) 

studied by Chapman belong to this category, although even 
these may have experienced shielding to some lesser degree. 

In this regard. Chapman's calculations may provide a lower bound 
on tektite entary into the atmosphere. On the other hand, the 
existence of specimen with barely ablated sturfaces precludes 
the possibility of entry at speeds much greater than 11 km/sec. 
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